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Abstract 

We establish the variational principle of Kolmogorov-Petrovsky- 
Piskunov (KPP) front speeds in temporally random shear flows inside 
an infinite cylinder, under suitable assumptions of the shear field. A 
key quantity in the variational principle is the almost sure Lyapunov 
exponent of a heat operator with random potential. The variational 
principle then allows us to bound and compute the front speeds. We 
show the linear and quadratic laws of speed enhancement as well as a 
resonance-like dependence of front speed on the temporal shear corre- 
lation length. To prove the variational principle, we use the compari- 
son principle of solutions, the path integral representation of solutions, 
and large deviation estimates of the associated stochastic flows. 



1 Introduction 

Reaction-diffusion front propagation in strongly time dependent random me- 
dia arises in premixed flame propagation problems ([13123 ED 1221231201 and 
references), interacting particle systems (|2S1III1 and references) and popu- 
lation biology (122] and references). A fundamental issue is to characterize, 
bound and compute the large time front speed, an upscaled quantity that 

*Department of Mathematics, University of Texas at Austin, Austin, TX 78712 
(jnolen@math.utexas.edu). 

^Department of Mathematics, University of California at Irvine, Irvine, CA 92697. 
(jxin@math.uci.edu). 



1 



depends on statistics of the random media in a highly nontrivial manner. 
In combustion hterature, ad hoc and formal procedures abound for approx- 
imation, such as closures and renormalization group methods EH] • In 
this paper, we establish variational principles of propagation speeds of KPP 
react ion- diffusion fronts through temporally random shear flows inside an 
infinite cylinder. The variational characterization then allows us to estimate 
and compute the statistical properties of front speeds with both accuracy 
and ease. 

The model equation is: 

Ut = ^A,u + B ■\/,u + fiu), (1.1) 

where u = u{z,t), z = {x,y) E R x Q, Q a, bounded open subset of K^~^ 
with Lipschitz continuous boundary, n > 2; B = t), 0, . . . , 0), b{y,t) 
is a stationary Gaussian process in t, with a given deterministic profile in 
y, to be made more precise later. The nonlinear function f{u) G C^([0, 1]) 
is a KPP nonlinearity: f{u) > for m G (0,1), /(O) = /(I) = 0, /'(O) = 
sup^g(o,i) fi'^)/'^- All example is f{u) = u{l — u). 

For compactly supported initial data bounded between and 1, solutions 
of (jLip develop into propagating fronts separating the cylindrical domain 
into a region where u ^ 1 and the rest where m ~ 0, which correspond to 
burned (hot) and unburned (cold) states in combustion. In case B is periodic 
in z and t, KPP type front dynamics and speeds have been recently studied 
for both shear and more general incompressible flows [201 1221 1211 12Z1 1211 120] ■ 
Exact traveling front solutions exist |2Z1 12H1 12E1 , extending those in spatially 
periodic media, [H |31 EHl iHj , see also jS] and [l?] for reviews. 

For temporally random shear flows, it is more efficient to study front so- 
lutions asymptotically without constructing exact traveling fronts. This line 
of work goes back to Freidlin [14jj where variational principles of KPP front 
speeds in spatially periodic media are obtained by combining the large devia- 
tion techniques and Feynman-Kac representation formulas of KPP solutions. 
We shall further develop this approach to treat the temporally random shear 
flows which generate more complexities in path integrals and unbounded 
variations in time. 

Let us make precise our assumptions on the shear field. The function 
b{y,t) = b{y,t,Lj) is a mean zero Gaussian random field over {y,t), periodic 
in y with period L for each fixed t, and stationary in t for each fixed y. The 
field b is defined over probability space {Cl, JF, Q) and has covariance function 
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T{yi,y2,ti,t2) = EQ[b{yi,tt)b{y2,t2)]. The following assumptions hold on 

Ky,ty. 

Al: (Periodicity in y) Let Cp^{D) denote the space of Lipschitz continuous 
functions that are periodic on the period cell D = [0,L]"~^. For each 
(1) E fl, there is a continuous map J{-,lj) : [0, +00) —>■ Cp^{D) such 
that b{-,t,(jj) = J(t,u). 

A2: (Stationarity in t) For each s G -R"*" there is a measure preserving trans- 
formation Tg : Cl ^ Q such that b{y, ■ + s,uj) = b{y, ■,TsLj). Hence, F 
depends only on yi, y2 and |t — s|. 

A3: (Ergodicity) The transformation Tg is ergodic: if a set A e JF is invari- 
ant under the transformation r^, then either Q{A) = or Q{A) = 1. 

A4: The field b is mean zero, almost surely continuous in {y,t), and has 
uniformly bounded variance: 

E[b{y,t)]=0 E[b{y,ty]<a^ for all y E D,t > 0. (1.2) 
A5: The function r(r) = sup^^ r(?/i, 1/2, 0, r) is integrable over [0, 00): 

POO 

/ T{r) dr = Pi < 00 (1.3) 
Jo 

for some finite constant pi > 0. 
A6: There is a finite constant P2 > such that 

\T{yi,y2,s,t) - F(?/i,?/3,s,t)| < - 2/2|f (|s - t|). 

For example, a field satisfying Assumptions A1-A6 might have the form 
b{y,t) = '^f^ib{{y)b{{t), where the functions b{(y) are deterministic, Lip- 
schitz continuous and periodic over D, and the functions &2('^) mean 
zero, stationary Gaussian fields in t. 

Before stating the main results, let us define the family of Markov pro- 
cesses associated with the linear part of the operator in (jl.ll) . For a fixed 
CO eil and for each 2 G i?", t > 0, let Z''\s) = {X''\s),Y''\s)) G i?" solve 
the Ito equation: 

dX'^\s) = dWi{s) +b{W2{s),t- s)ds 

dY''\s) = dW2{s) (1.4) 
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with initial condition Z^'*(0) = z = {x,y) e i?", where W{s) = (W^i(s), ^^(s)) G 
i?" is the n-dimensional Wiener process with W{0) =0. Let P^'* denote the 
corresponding family of measures on C([0, t]; i?"). As we will see, the KPP 
front speed depends on large deviations of the random variable 

Vlint) = "-^^^^ (1.5) 

which is the first component of the average speed of a trajectory over time 
interval [0, nt]. We consider only the first component since we are concerned 
only with propagation in the x direction. The need for the parameter k 
results from the time dependence of the field b{y, t) and will become more 
apparent later. 

Now we state the main results. First, the following lemma allows us to 
characterize the speed of propagation: 

Lemma 1.1 Assume that A1-A6 hold for the process b{y,t). Then for each 
A G -R, the limit 

/i(A, z) = /i(A) = /'(O) + lim - log E [e-^^^''^] (1.6) 

exists and is a finite constant, almost surely with respect to measure Q and 
independent of z E -R". Moreover, /i(A) is a convex, positive, and even 
function of X. Also, yu(A)/|A| — > +oo as |A| — > oo. 

If we let H{c) be the Legendre transform of /i(A) 

H{c) = sup[c • A — /i(A)], 
Aei? 

then we find that the speed of propagation can be bounded above in terms 
of if. 

Theorem 1.1 (Upper bound on front speed) Let b{y , t , uj) satisfy assump- 
tions Al - A6. Let u{x, y, t, lj) solve with initial condition u{x, y, 0, u) = 
Uo{x), where Uo{x) G [0,1] has compact support and is independent of u. 
Then, for any closed set F G {c E R\ H{c) > 0}, 

lim sup u{ct, y, t,uj) = 
uniformly in c E F, for almost every a; G ^2. 
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Therefore, if we define the constant c* > by the variational formula 



c' = M^, (17) 

A>0 A 

we see from the definition of H that the front spreads asymptotically with 
speed no greater than c* in the positive x direction and with speed no greater 
than c* in the negative x direction. Although the solution u depends onuj & Q 
since 5 is a random variable over Cl, the function H{c) and the speed c* are 
independent of a). They are almost surely constant with respect to Q, a 
consequence of the ergodicity assumption A3. We will generally suppress the 
dependence of m on for clarity of notation. 

The constant c* is also almost surely a lower bound on the speed of 
propagation: 

Theorem 1.2 (Lower bound on front speed) Let b{y,t,Lj) satisfy assump- 
tions Al - A6. Let u{x, y, t, uj) solve with initial condition u{x, y, 0, lj) = 
Uo{x), where Uo{x) G [0,1] has compact support and is independent of uj. 
Then, for any compact set K d {c & R\ H{c) < 0}, 

lim inf u{ct,y,t,uj) = 1 

t~*oo yGD 

uniformly in c & K , for almost every uj eCl. 

Theorems 11.11 and 11.21 extend our recent results on KPP front speeds 
in temporally periodic incompressible flows ^28j i26j and a classical result 
of Freidlin (Theorem 7.3.1, p. 494 of J3]) where he treated the case of 
spatially periodic advecting flows. Our proofs are built on his, with additional 
ingredients to handle both the time-dependence and the stochastic nature of 
the field B. For example, in the periodic case, /i(A) is the principal eigenvalue 
of a periodic-parabolic operator 28j26,, and perturbation theory |17j implies 
that /i(A) is differentiable in A. It then follows from Theorem 7.1.1 and 
Theorem 7.1.2 of jT3] that the random variables ri\{t) satisfy a large deviation 
principle with convex rate function S{c) given by (j4.3|) . However, if /i(A) is 
not known to be differentiable, the large deviation property needs a new 
proof. In the present case, /i(A) is not an eigenvalue of a linear operator, so 
we cannot readily apply the perturbation theory 17J to get differentiability. 
Instead, we will show that a rate function exists and is convex, and that it 
continues to satisfy (j4.3p . In fact, /i(A) is related to the almost sure principal 
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Lyapunov exponent of a heat operator with random potential (HS] , known as 
the parabohc Anderson problem (fOl^] and references). Dynamical aspects 
of principal Lyapunov exponents as an extension of principal eigenvalues are 
recently studied in [21]. Regularity of /i(A) is an interesting problem in itself. 

The paper is organized as follows. In Section 2, we prove Theorem 11.11 
In Section 3, we adapt Freidlin's method to prove Theorem 11.21 assuming 
a few technical estimates as well as a crucial lower bound stated in Lemma 
13.11 In Sections 4 and 5 we will prove the lower bound of Lemma 13.11 and 
necessary large deviation estimates. It is in these sections that we resolve 
the main difficulties that result from the additional stochastic time depen- 
dence. We will make frequent use of the subadditive ergodic theorem and 
the Borell inequality for Gaussian fields [Ill2llin]- In Section 6, we use the 
variational formula ()1.7p to numerically compute the front speed c* in an 
example where the temporal random process is Ornstein-Uhlenbeck (0-U). 
We study the dependence of front speeds on the covariance of the random 
media as well as their growth laws in small and large advection limits. The 
speed grows linearly in the limit of large advection, and quadratically in 
the limit of small advection. The variational principle also yields analytical 
bounds on c*, demonstrating the enhancement by random shear fiows and 
that small temporal correlation reduces the front speed. This is analogous 
to the decrease of front speeds with increasing frequency of temporally os- 
cillating periodic shear flows 1771 12H] • Furthermore, there is an optimal 
correlation length that maximizes the enhancement. Linear and quadratic 
speed growth laws are known for deterministic flow patterns with channel 
structures (0 IH El UHl UH El] and references), also for spatially random 
shears inside infinite cylinders [211 (SO] or white in time Gaussian shears in the 
entire space [SH]- The speed growth laws known to date are not sensitive in 
the form of nonlinearities as long as fronts propagate out of the initial data. 
In this sense, KPP plays the role of solvable model and KPP front speeds 
carry universal properties. In Section 7, we make concluding remarks. Fi- 
nally in the two Appendices, we prove Lemma f 1.1 1 and several key technical 
bounds. 

Though the arguments in our proofs rely on the periodicity of t) in 
y to provide compactness in the y dimensions, they can be easily modified 
to solve the same problem in an infinite cylinder with the zero Neumann 
boundary condition on the sides of the cylinder. The compactness property 
remains, and the process X^'*(s) just needs to be refiected when it hits the 
boundary R x dVt. It is also not necessary for the process h{y, t) to be 
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Gaussian. Our proofs of the bounds in Sections 4 and 5 shall rely on the 
powerful Borell inequality for Gaussian process, yet it is easy to see that if 
the estimates of Lemma 13.21 and the Appendices hold for a given process, 
and if assumptions Al - A6 hold, the main results extend. 



2 Proof of Theorem 11.11 

The proof of Theorem 11.11 is based on the assumption that f{u) < f'{0)u. 
This allows us to construct a super-solution to equation ()1.1|1 as follows. Note 
that the integral form of ()1.4j) is 



X + [ b{W^{T),t-T)dT + W^{s) 

Jo 



(2.1) 



where W^{s) denotes a Wiener process starting at W^{0) = z, P-a.s. Using 
()2.H) . we can express ()1.6|) as 



/'(0)+ lim-logE, 



/'(O) + 



lim - log Ey 

t—*oo t 



(2.2) 



Here we have used the independence of Wi and W2- Note that we will 
sometimes use p(A) = /i(A) — /'(O) — A^ to refer to the limit on the right hand 
side of (|2.2p . Using the Feynman-Kac formula, we see that yu(A) is just 



lim - log Ey 



X\^/2+f'(Q)~\b(Wy{s),t~s)ds 



lim -log$(?/,t) 



where $ = $(?/, t) > is periodic in y and solves the auxiliary initial value 
problem 

= iA,<|. + (AV2 + /'(O) - A6(y, 

$(y, 0)^1. 

Suppose c > c* (i.e. H{c) > 0). Then for e > sufficiently small, there 
exists A > such that yu(A) < (c — e)A. By Lemma fl.ll /ifA) is finite, so that 



lim -log<l>(?/,t) =/i(A), 



and there is a function R = R{y,t) such that |-R| ^ as t ^ oo, uniformly 
in y, such that 

Defining ip^x.y.t) = e^^^^{y,t), we see that ip solves the equation 

^(a;,y,0) = e-^^ 
By the properties of /, we then have 

when ip G [0, 1]. Multiplying ip hj a constant C if necessary, we may assume 
ip{z, 0) > ^0(2;), and by the maximum principle we must have u{z, t) < ip{z, t) 
for alH > 0. If we define the half-space 

S+ = {{x,y) e i?" I X > r} 

we now see that 

lim sup u{x,t) < lim sup ip{z,t) 



= lim sup e (f){y,t) 

< lim sup e-^=*e^(^~^)*+^(^'*)* = 

since \R{y,t)\ < y for t sufficiently large. A similar argument holds for 
propagation in the —x direction. This completes the proof of Theorem 11.11 
□ 



3 Proof of Theorem HH 

Proving Theorem 11.21 is more difficult, and we will need the following impor- 
tant estimates. The first is a lower bound analogous to Lemma 7.3.3 of P^ : 
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Lemma 3.1 For any compact set K G {c E R\ H{c) > 0}, 

liminf-log inf u{ct,y,t) > —ma:xH{c). (3-3) 

t^oo t ceK,yeD ceK 

The second estimate gives a coarse bound on very large excursions of the 
random process X^'*: 

Lemma 3.2 There are constants Ki,K2 > independent of k, E (0, 1] such 
that, except on a set of Q -measure zero, 



sup P sup -x\>r]t \ < Kie-^^''^''^ 

z<^R" \ se[o,Kt] j 

for any k G (0, 1], ?7 > 0, and for t sufficiently large depending on uj, k, and 
V- 

These lemmas represent the main technical difficulty in extending the 
work of to the present case with a stochastic time dependence in b{y,t). 
For the moment, however, we delay the proof of these lemmas and show how 
these bounds lead to Theorem 11.21 Lemma f3. II is proved in the next section; 
Lemma 13.21 is proved in the appendix. The proof of Theorem 11.21 is based 
on the observation that when u < h < 1, the reaction rate can be bounded 
below. For each u G (0, 1], define the reaction rate ( by 

C(.) 



u 

and C(0) = /'(O). Now equation ()1.H) can be written 

Ut = ^A^u + b{y,t)u^ + C{u)u. (3.4) 

By the properties of f{u) we see that C{u) > for m G [0,1), ({u) is contin- 
uous for u G [0, 1], and C(0) > C(''^) ^oi any u G [0, 1]. If /i G (0, 1) we define 
a lower bound on (: 

a= inf CH>0. 

mG{0,/i) 

So, in regions where u is bounded away from one, the reaction rate can be 
bounded below by C/i > 0. 
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For a fixed a) G f2, we can estimate u{z, t) using the Feynman-Kac formula 
for the solution of (j3.4|) : 

u{z, t) = E [e-^o at-s,u{z-^'\s),t-s))ds^^^^z,t^^^^A ^ ^3_5^) 

where the expectation is with respect to measure P^'*. If r is any stopping 
time, we also have 

u{z,t) = E L/o^^f(*-^'-(^^'*W.t--))'^-M(Z^.*(t - (t A r)),t - (t A r))j , (3.6) 

where t/\T = min(t, r). Therefore, we can obtain estimates on u by carefully 
choosing stopping times and restricting the expectation to certain sets of 
paths. The exponential term inside the expectation will be large when the 
path passes through regions where u is small and the reaction rate is 

large. On the other hand, if u{Z^'^(t — {t A T)),t — (t A r)) is too small, then 
the expectation whole may be small. 

Now we present the argument of (see p. 494). For s G i?, define the 

set 

*(s) = {ceR\ H{c) = s} and ^{s) = {c E R\ H{c) < s} 
and then for any 5 > and T > 1, 

[{l}x^(5)]U[ U {t}x{m{5))]\xR-~\ 

, l<t<T J 

This defines the boundary of a region that spreads outward in x, linearly in 
t. Outside this region u may be close to zero, but on the boundary of this 
region, we have the crucial lower bound from Lemma f3. II 

m(x, y, s) > e~2^* for all (x, y, s) G Tt (3.7) 

for t sufficiently large. 

For h G (0, 1), and z fixed, t, r/ > 0, define the Markov times 

ahit) = min{s G [0,t]| M(Z^'*(s),t - s) > /i} 
arit) = mm{s e[0,t]\{Z''\s),t-s) eVt} 
T^it) = mm{s G [0,t]\ |X^'*(s) - x\ > T]t} 

(We set these variables equal to +oo if the sets on the right are empty.) Then 
P^'^{o'h{t) < t) is the probability that a particle starting at z will encounter 
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the "hot region", u > h, at or before time t. Since C ^ O5 it is clear from 
(jSSD that 

uict,y,t) > hP'''y'\crhit)<t). (3.8) 

For simphcity, we will write P*^*'* or just P'^^ to denote p('=*'2^)'* when the y 
and if: dependence is clear. If we choose c too large, we expect that this 
probability will be small, since x = ct would be far beyond the spreading 
front, beyond the region where u > h. On the other hand, we want to 
show that the probability is large foic&K (i.e. |c| < c*), for if we can 
show that for each h, P^^'^{ah{t) <t)^l uniformly over c E K, then ()3.8p 
implies the desired result: u{ct,y,t) —>■ 1 uniformly over c E K. Note that 
P'''\ah{t) >t) = P''''\ahit) > t, r,(t) > t) + P''^\ah{t) > t, r^{t) < t). By 
Lemma 13.21 

sup P^'* (r^(t) < t) ^ (3.9) 

as t 00 except on a set of Q— measure zero. So, it suffices to show that 

P"*'*((Tft,(t) > t, Triit) >t) ^0 (3.10) 

uniformly in c G -ft'. Note that in [T^, estimate ()3.9p followed from the 
uniform boundedness of the field B, a property that we do not have in the 
stochastic case. The rest now follows exactly as in ^3]. Choosing rj small 
and then r G (0, 1) sufficiently small, we have 

rt < arit) <t-l 

whenever r^(t) > t. Now, 

P'^'^'Mt) > t, T,{t) >t)< P^*'*(A) 

where A is the set A = {u G Q\ rt < ar < ah{t)}. This probability is 
bounded above by 



P"*'*(A) < 



ct 



Here we have used the crucial lower bound ()3.7|1 . The expression on the right 
is bounded by 



ct 



j;^ Ciu{z,t-s))ds^-irta^^^2{crr),t - crr)Y/\A 



11 



Then using Holder's inequality, we have 

P^*'*(a;,(t) > t, T^{t) >t) < e^*[M(ct,t)]^/2g-iChrt 

< e*(^-H'>-) (3.11) 

The last term goes to zero if S is sufficiently small, depending on h. This 
completes the proof of Theorem I1.21 □ 

Note that the only difference between this argument and that of Freidlin 
in the periodic case is the manner in which ()3.7j) and ()3.9|1 are obtained. 



4 Proof of Lemma 13.11 

The main issue in proving the estimate of Lemma 13.11 (and thus the lower 
bound ()3.7|l ) is whether the random variable 

Vli-t) = "-^^^ (4.1) 

satisfies a large deviation principle with a convex rate function that can be 
characterized by /i(A), almost surely with respect to Q. The variable rj\{nt) 
is the first component of the average speed of a trajectory over time interval 
[0, Kt]. 

Definition 4.1 For fixed uj G VL, the random variables riKnt) satisfy a large 
deviation principle with a convex rate function S{c) if there exists a convex 
function S{c), independent of z E i?", such that 

(i) For each s > 0, the set $(s) = {c G -R| S{c) < s} is compact. 

(a) For any 6,h > 0, there exists t^ > such that for all t > t^ 

P{d{r^l{Kt),^s))>5)<e-'^^-^\ 

(Hi) For any 5,h > 0, there exists t^ > such that for all t > t^ 

P {r]l{Kt) G Us{c)) > e-'^*(^W+'*). (4.2) 

If such a function S{c) exists, it might depend on the parameter k, G (0, 1], 
and it might depend on (Jj E Cl. However, we will show that 
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Proposition 4.1 Suppose that assumptions A1-A6 hold. Then almost surely 
with respect to Q, the random variables f]l{Kt) satisfy a large deviation prin- 
ciple (with respect to P) with a convex rate function S{c) that is independent 
of K, and a; G . 

We postpone the proof for the moment while we finish the proof of Lemma 
13.11 By Lemma fl .11 the quantity 

is well defined and is almost surely constant with respect to Q for A G -R, 
independently of k. Since, by our assumption of Proposition l4.1l the variables 
f]l{Kt) have a convex rate function, it follows that S{c) has the following 
characterization (see Section 5.1 of jT3]): 

5(c) = sup[cA -^(A)]. (4.3) 

Ae-R 

This characterization does not hold if S{c) is not convex. Let us emphasize 
that S{c) is independent of k G (0, 1] and lj E Q, although the constants to 
in Definition 14. II mav depend on k,,lj. 
Now, by definition of S{c), 

liminf — log inf P WAi^t) G Us{c)] > -S{c) > -oo, (4.4) 

and the lower bound ()3.3j) of Lemma 13.11 can be written 

liminf-log inf u{ct,y,t) > f'{0) — max S{c). (4.5) 

t — >oo t c<^K,yeD ceK 

To prove the lower bound we now use the Feynman-Kac formula to relate 
fl4.5j) to ()4.4j) . as in the arguments of Freidlin in Lemma 7.3.2 in jT3]. The 
compactness of K implies that it suffices to show that given any e > 0, and 
any c for which H{c) > 0, 

liminf (-log inf u{ct,y,t)] > f'{0) - S{c) - e (4.6) 

t^oo yt c€Us{c),y£D J 

for 5 > sufficiently small. Without loss of generality, we assume that 

Uo{x) > XUs{o){x) (4.7) 
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for S sufficiently small. That is, Uq = 1 whenever, |x| < S. We define q to be 
the limit on the left hand side of (j4.6|) : 

g = lim inf ( - log inf u{x,y,t) 

t->oo yt x€Ust{ct),y£D 

Without loss of generahty, we assume < c* < c where c* is defined by fjl.7|) . 

Step 1: The first step is essentially the same as in jT3]. Suppose for the 
moment that we known q is finite. By the representation (|3.6p we have for 
any n G (0, 1] 

inf u{tc, y, t) > inf E \ed(i-'Mz(s),t-s))ds^^2^^ -Kt),t- Kt)xA 

c&Us{c),y&D c&Usic),y&D L 

(4.8) 

for any set ^s<rmeasurable set A. Recall that when u < h, the reaction 
rate ({u) is bounded below by C/i > 0. If we choose A to be the set of paths 
satisfying both 

X^'*(fi;t)Gf/(i_,)5t((l-K)te) (4.9) 

and 

u{Z'^\s),t- s) < h for all s G [0,Kt], (4.10) 
then from ()4.8|) and the assumption that q is finite we have a lower bound 

q > Cfe + liminf — log inf P(A), (4.11) 

t^oo Kt ceUs(c),y€D 

provided that the limit on the right also exists and is finite. 

Step 2: Now we bound the right hand side of ()4.11|) and show how 
it relates to ()4.4|) . Let 6 be sufficiently small so that we can pick c' with 
c* < c' < c — 66. By Theorem 11.11 for any h G (0, 1) there is a constant 
to > such that 

u{x,y,t) < h for all x > c't, y G t > tg- 

Now if K < 1/2 and 

sup \X^'\s) -{t- s)c\ < 3St, (4.12) 

then X^'*(s) > c'(t — s) for all s G [0, Kt]. Thus, fl4.10|) is achieved along such 
paths when t > 2tQ. Next, if c G Us{c) is written c = c + 6ei for some ei with 
|ei| < 1, then define c = c + 26ei. Then for any |e2| < 1 

ct — K.tc + Kt5e2 G ?7(i-K)5t((l — K)ct). (4-13) 
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It follows that for each c G Us{c) there is a c G f/2<5(c) such that (|4.9p is 
achieved whenever ri\{Kt) G Us{c), where rj is defined by (|4.ip . This gives 
us a lower bound on P{A) in terms of the riKut), the average speed of a 
trajectory over [0, nt]: 

inf P{A) > (4.14) 

ceUs(c),y&D,z=ct 

inf P I sup |X"'*(s) - (t - s)c| < 35t, r]\{Kt) G ^/^(c) 

For K sufficiently small, k, < (25)/(3max(l, |c|)), we see that 

sup P I sup \X''\s) - (t - s)c| >36t] < sup P ( sup \X''\s) - x| > (5t/3 

By Lemma 13.21 there are constants Ki,K2 > independent of k such that 
(except possibly on a set of Q-measure zero) 

sup P I sup \X''\s) -x\> St/s] < Kie-^^^^''^ (4.15) 

for t sufficiently large depending on Cj. Therefore, for any M > 0, by choosing 
K arbitrarily small, we can make K^S^ j > M, so that 

lim sup — log( sup P ( sup |X'''*(s) - (t - s)c\ >26t])< 



t—*oo 



< lim — log(irie-^2<52i/«) < _M. 
Therefore, from ()4.14j) and ()4.1H) we now see that for k sufficiently small, 

q > Ch + liminf— inf P (r]iUt) e Us(c)) (4.16) 

t^oo nt ceU2s(c),zeR" 

provided that the limit on the right is finite and bounded below, indepen- 
dently of K. However, this follows immediately from Proposition 14.11 and the 
lower bound (jOj)- Then (gH) follows by letting h ^ so that (h ^ /'(O). 

Step 3: It remains to establish the initial claim that q > — oo. To see 
this, define for any Ci G P 

q5{ci,t)= inf P4X^'*(t) G f/5(0)) , (4.17) 

x&U6itci),y(^D 
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a random variable over Let use write X^'^'*(s) as 

X''\s) = X + P'\s) + Wi{s) 

where /^'* is the first integral term in ()2.H1 and Vri(O) = 0. For uj ^ Q fixed, 
let M > and define the set 

Am = {w\ sup \P{s)\<Mt}. 

yeD,se[o,t] 

Using the fact that Wi and W2 are independent, we derive the lower bound 
for X G Us{cit) 

P{Xt' G Usm > P{W,it) G Us{0) - Pit) - x; A) 

> PiWiit) G UsiO) + Mt + \ci\t + 6)PiAM) 

> ^ P{Am). (4.18) 

By Lemma ini2l -P(^m) > 1/2 for t sufficiently large, depending on a) and 
M. Therefore, there is a finite constant Ki > depending only on ci and M 
such that ^ 

liminf -logg5(ci,t) = -i^i (4.19) 

t^oo t 

is finite almost surely with respect to Q. Now for ci < C2, we get two 
constants Ki and ^^^2 such that the limit (j4.19p holds almost surely with 
respect to Q. Let Ci = c — 6 and C2 = c + 6 and define the Markov time 

n{z,t) = inf {s > 0| X'^'^s) G (-cx),Ci(t - s)] U [c2(t - s), 00); or s > t| 

which is the first exit time of the process X^'*(s) from the region defined by 
{{x,y,T) \ X G {ciT,C2t),t > 0,y G D}. By formula we have 

inf u{ct,y,t) = inf u{x,y,t) (4.20) 

ceUs{c),yeD xe{c-it,C2t),yeD 

> inf E [M(Z^'*(7r),t-7r 

2:e(cit,C2t),ye-D 



> inf m(ciS, y, s) A inf u{c2S,y,s] 

\s€[0,t],y(^D J \s(^[0,t],yeD 

It follows from ()4.7|) and ()4.19p . that there are constants Ci, C2 > (depend- 
ing on u) such that 



u 



CiS,y,s)>Cie'^^'' and u(c2S, s) > Cse-^^^* (4.21) 
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for all s > 0. Combining (IOn|l with (jOT|l . it is clear that q > -2(7^1 Vi^z) > 
— oo. Having shown that q is finite, we have completed the proof of Lemma 

o □ 



For later use, let us show that for all t > 0, log(g5(c, t)) is integrable with 
respect to Q. Using (|4.18|) we see that 

1 1/(5 (Al"t+|ci |t+25)2 \ 

-logg,(c,t) > -log(^-^e ^ P{AM)j (4.22) 

> -C^ + -logie ^ P{Am: 

for a constant Ci > independent of c, for t > 1. Let g be the term inside 
the logarithm: 

(Mt+|ci |t+2i)2 

9 = e i P(Am). (4.23) 

Then for a > 2Ci, 

Q (^hogqs{c,t) < -a^ < g01og^<-a/2 

= Q (p(Am) < e-^/^e^^^^^^^) .(4.24) 

Also, from Lemma f9.2t 

Q (PiAM) < 1 - e"^^*^'*/2j < i^^e-^^*^'*/^ (4.25) 

With a little algebra, one can see that there exist constant K^, > 
independent of t such that whenever t > 1, M = K'j^^fa^ and a > K^c^ , we 
have 

^-^m^iMlM^ < 1/2 < 1 - e-^^*^^*/l 
By combining ()4.24|1 and ()4.25p . we now conclude that 



Q Qlogg5(c,t) < -a^ < i^ie-^^^'^^/^ 



(4.26) 
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whenever a > Ki(? and t > 1. It follows that for t sufficiently large, inde- 
pendent of CJ, 



[I ^ log <?5(c,t) 



= j Q (^j\ogq5{c,t)\> da 

POO 

< Ki^+ / irie-^^^3"*/2 < oo. (4.27) 



This is bounded uniformly in t, for t > 1. 

5 Proof of Large Deviation Estimates 

In this section we prove Proposition 14.11 We work first with the case k = 1. 
Because b is independent of x, it suffices to show that the proposition holds 
for X = 0. Define for c & R and r < s <t 

qlic,r,t)=Pi-X''y'\t-r) E Us(t-r){c{t - r))) = P{r,l{t - r) E Us{c)) 

and 

(lsic,r,t) = supgf(c,r,t) 
y 

(h{c,r,t) = inf g|(c,r,t). 

y 

We will use the subadditive ergodic theorem to show that (1/t) logg^(c, 0, t) 
and (l/t) logg/(c, 0, t) converge almost surely to a finite constant. Define 
the events 

A = {-X'^y^^t - r) E Usit-rMt - r))} = {T^lit - r) E Us{c)} 
B = {-X''y'\t-s)EUs^t-sMt-s))} = {vl{t-s)EUs{c)}. 

Note that event B is J^s<t measurable for t >t — s, where J^r is the cr-algebra 
generated by (Z^'*(s))s<t-. By the Markov property of the Wiener process we 
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y 

inf 

y 

inf 

y 

inf 

y 

> infE 

y 



have: 

qj{c,0,t) = MP{A)>miP{AnB) 

y y 

> inf P({-X°'^'*(t - r) + X^'^'^t -s)e f/5(r-s)(c(s - r))} n B) 
= infE[xBP[{-X°'^'*(t-r)+X°'^'*(t-s)Gf/5(,„,)(c(s-r))} \ J^t- 

inf E [xBP[{-X°'^'*(t - r) + - s) G f/5(.„.)(c(5 - 0)} I " s 

Xb MP[{-X^'''%s -r)e Usis^rMs - r))}] 
= inf P({-X°'^'^(s - r) G f/5(.-.)(c(s - r))}) inf P(5) 
= (lj{c,r,s)qj{c,s,t). 

Therefore, log(g^ (c, s, t)) is super-additive for each c E R. By the stationarity 
of b, 

Thq;{c,r,t) = Tf,mfP{-X^'y'\t-r)eUsit-rMt-r))) 
y 

= infPi-X'^y^'^^t - r) G Us^t-rMt - r))) 
= qj{c,r + h,t + h). 

For any e > 0, we can bound q below by translating in x and using (j4.17j) : 

qj{c,r,t) > Trq,{c,t- r) = inf P(X^'^'*(t - r) G t/e(cr)) (5.28) 

ye-D,a::el/E(ct) 

for t — r sufficiently large. Hence, log(g^(c, r, t)) is integrable by ()4.27|) . 
Kingman's ergodic theorem JH] now implies that the limit 

1 _ 1 _ 

lim -logg^ (c,0,n) = sup -Pgflogg^ {c,0,n)] = -Ss{c) (5.29) 

n^oo 77, „>o n 

exists and is a finite constant, Q-a.s, because of the ergodicity assumption A3. 
Using an idea in |2], we can now extend convergence in (j5.29p to continuous 
time. Let 

g{uj) = sup \\og{qj{c,r,t))\. 

r,te[0,2] 
\r~t\>l 
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Let T(cj) = supyg£i_4g[Q 2]|&(y, Then for all < r < t < 2, 

pt—r 

sup| / h{Wy{T),t-T)dT\<T\t-r\. 
yen Jo 

As in dHHl), 

P{-X''y'\t - r) G ?75(t_,)(c(t - r))) > P{Wiit - r) e Usit-r){0) + T(t - r) + |c|(t - r)) 

6\t — r\ {t-rf(r+\c\+S)'^ 
> ' ' e 2(t-r) 

^27r|t - r| 

Therefore, 

for some constants Ki,K2 > that depend on 6 and c. Hence (7((u)) is 
integrable with respect to Q. By the super-additivity of qg{c,r,t), 

q^{c,0,n- 1) - r„_i5( < qj{c,0,t) < qj{c,0,n + 2) + r^g (5.30) 

whenever t & [n, n + 1) , n ^ Z . As in the proof of Theorem 2.5 of [2], one 
can show that ^t„(7 — > almost surely as — > oo since (7 is integrable. It 
now follows from ()5.30|) that 

lim ylogg7(c,0,t) = -Ss{c) 

almost surely with respect to Q. 

For each c E R, Ss{c) can be bounded above independently of 5 > using 
fl5.28|) and ()4.19|) . From the definition, it is clear that Ss{c) > for all 6, and 
Ssi{c) < Ss^ic) whenever 61 > 62. Therefore, we define for each cE R 

S{c) = lim Ss{c) = sup 5*5(0) G [0, +00). 

Lemma 5.1 For all 6 > 0, the functions Ss{c) are continuous and convex 
in c. Also, S{c) is continuous and convex in c. 

Proof: The continuity and convexity of S{c) follows immediately from the 
continuity and convexity of Ss{c), since the functions Ss{c) converge pointwise 
to the finite limit S{c). The convexity of Ss{c) follows from the Markov 
property of the process X^'*, as follows. 
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Let p G (0, 1) and Cq = pci + (1 — p)c2- Let t > and denote ti = pt, 
t2 = (1 — p)t. Then we see that 

g,-(co,0,t) = infP(-X°'^'*(t)Gf/5t(cot)) 

> inf P (-XO'^'*(t) G f/5t(cot) , -XO'^'*(ti) G f/5t,(citi)) 

?/ 

> inf P (-X0'^'*-*Ht2) e f/5t,(c2t2)) inf P G UstMh)) 

y y 

= (Ij{c2,0,t2)qj{ci,t2,t). 

Hence 

-j^ogqj{co,0,t) < -^Iogg7(c2,0, (1 - ^logg7(ci, (1 -p)t,tX5.31) 



^logg^ (c2,0, (1 - r(i_p)t Qlogg^ (ci,0,pt)^ 



By the stationarity of b, the random variable 



converges in distribution to pSs{ci). Therefore, there is a set M C with 
Q{M) = such that for each u & Cl\M and each e > 0, there is an increasing 
sequence tj — > cxo as j ^ oo, such that for j sufficiently large, 

-ni-p)t, (^^logg7(ci,0,ptj)^ <pSs{ci) + e. (5.32) 

Note that the other two terms in ()5.31|) are random variables that converge 
to constants (Q-a.s.) as t oo. For fixed u E Q and e < 0, we can pick t 
sufficiently large so that 

-^logg,(c2,0,(l-p)t) < il-p)Ssic2) + e, (5.33) 

and 

-^logg,-(co,0,t) > Ss{co)~e. (5.34) 
Now by (fOT|l and ((OS)- (ESI we have 

Ssico) - e < (1 - p)Ss{c2) + pSs{c^) + 2e. 
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Since e > was arbitrarily chosen, we infer that 

Ssico) < {l-p)Ss{c2)+pSs{c,). (5.35) 

So, Ss{c) is convex and must also be continuous in c, since it is finite for 
every c G -R. □ 

This establishes the existence and convexity of the function S{c). Part 
(iii) of the Definition 14 . II follows from the definition of Ss{c) and the fact that 
Ssic)/'Sic). 

To finish the proof of Proposition 14.11 for k = 1, we must establish a 
Harnack-type inequality to relate the quantities 

P(-XO'^'*(t) G Ustict)) and P(-X°'J^''*(t) G Ust{ct)) 

for y,y' G D. This will allow us to remove the inf^ in the definition of q 
and S{c) and to establish parts (i) and (ii) of Definition 14.11 We prove the 
following lemma 

Lemma 5.2 There are constants Ki, K2, K3 > such that for all /t G (0, 1], 
c e R, e > 0, and 6 > 0, 

miP{r]l{Kt) G t/(i+.)5(c)) > KssnpP{vl{Kt) G Us{c)) - K^e-""''"'"^"'" 

^ z 

Proof: For clarity we let k = 1. Extension to /t < 1 is straightforward, as 
in the proof of Lemma (3. 21 in the appendix. First, note that 

= xo + + Wi{t), = 

and 

- = iy+w^(')^-'(t -s) + Wi{t) - Wi{s). 

With out loss of generality, we assume xq = 0. Therefore, 

p{x''y'\t) G Ust{-ct)) 

= PiX'^'y^t) - G Ust{-ct - 

= E [Pil^'^it -s) + Wi{t) - Wi{s) G Ust{-ct - a)|X°'^'*(s) = a, W2{s) = (3 
= E (P{l^^\t -s) + W,{t - s) G Ust{-ct - a)\X^'y'\s) = a, W2{s) = f3 - y) 
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Here we used Wi{t) to denote a Wiener process, Vri(O) = 0, P-a.s. Now we 
see that 

P{X'''y'\t) G Ust{-ct), |X°'^'*(s)| < M) 
< E (p(/^'*-^(t -s) + W,{t-s)e Ust+M{~-ct)\ W2{s) = /?-!/)) 

and 

p{x''y'\s) G Usti-ct), \x''y\s)\ < M) 

> E [P{l^^'-'{t -s) + W^{t-s)e Ust-M{-ct)\ W2{s) = P-y)). 

For any y' G -D, we have for fixed s G (0, t] 

P(XO'^''*(t) G Ust+2M{-ct), \X'^y''\s)\ < M) 

> E (p(/^'*-^(t -s) + W^{t - s) G Ust+M{-ct)\ W2{s) = /3 - y')) 

= / Ps{P-y')PiI''''-%t-s) + W^{t-s)eUst+M{~ct))df3 
Jd 

> Ci / ps{(3 - y)P{I^''-\t -s) + W^{t-s)e Ust+M{-ct)) d(3 

Jd 

= C^E [P{l^^'-'{t -s) + Writ -s)e Ust+M{-ct)\ W^is) = {3 - y)) 

>C,P{X'>'y^\t)e Usti-ct), \X^'y'\s)\ < M) (5.36) 

where Ps{r) denotes density function for distribution of W2{s) on the torus, 
and Ci is a positive constant depending on s. Using this inequahty, we see 
that 

p{x'''y''\t)eUst+2M{-ct)) = p{x''y'\t)eUst+2M{-ct), <M) 

+ p{x''y''\t) G f/,i+2M(-ct), \x'^y'\s)\ > M) 

> CiP(X°''^'*(t) G Usti-ct)) 

- c,p{x''y^'{t) G Usti-ct), \x'''y^\s)\ > m) 
+ p{x''y''\t) G Ust^2M{-ct), \x'^y''\s)\ > M). 

(5.37) 

For e > 0, let s = 1 and M = e5t/2 (M = e5/tt/2 when k < 1). It follows 
from the Borell inequality, as in proof of Lemma 13.21 that 

P f sup |X°'^'*(r)| > M I < K^e-""^"'" = Kie"^^^'^'*'/^ (5.38) 
\Te[o,i\,yeD J 
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for t sufficiently large. Now the lemma follows from (|5.37|) . □ 

Since e > is arbitrary, the lemma implies that 
Corollary 5.1 

Ss{c) = lim y logg7(c,0,t) = lim y logg/(c, 0,t) (5.39) 
almost surely with respect to Q. 

Now, using this estimate, we can establish parts (i) and (ii) from Defini- 
tion [O] From Lemma 13.21 there are constants Ki,K2 > such that for t 
sufficiently large, 

P{\vlit)\>c)<K^e- (5.40) 

This implies that lim|c|^oo 'S'(c)/|c| = +oo. Hence, $(s) is a bounded set, for 
each s > 0. Then by continuity of S{c), $(s) is compact. Let A be the set 

A = {ceR\ rf(r/*(t),$(s)) >(5}. 

We must show that for any fixed 6 > 0, h > 0, there is a to > such that for 
t>to, 

P{vl{t) eA)< e-'^'~''\ (5.41) 

Because of the bound ()5.4()|1 . it suffices to show that ()5.4H) holds with A re- 
placed by any compact subset A' of A (because K2C^ > s when c is sufficiently 
large). Pick e > small enough (at most e < 5) so that 

inf SJc') > s-h. (5.42) 

c'eA' 

This is possible since A' is compact and Ss{c') y S{c') as 5 — > 0, for any 
d E A'. Cover A' with a finite number of e-balls. Because the number of 
balls is finite, we conclude from ()5.39|) and ()5.42|) that for t sufficiently large 



P{r] e A') < e 



-t{s-h) 



This establishes Proposition 14. II in the case that k = 1. Now we extend the 
result to K < 1, as well. If /t G (0, 1) and 6 > 0, the stationarity of b{y,t) 
implies that for any z G -R" 

-^logP {4{nt) G Usic)) = -llogg|((l - K)t,t) ^ Ssic) (5.43) 

24 



in distribution (with respect to Q) as t — ^ oo, but this does not imply point- 
wise convergence. 

By definition, q^{s,t) < q^{s,t) < for all y,s,t. Also, as we have 

already shown, the Markov property implies sub and superadditivity: 

log qj (0, t) > log q^ (0, r) + log qj (r, t) 

and 

logg^(0,t) < logg^(0,r) +logg/(r,t) 
for all r, t > 0. In the same way, we also have 

logg,^(0,t) > log (0, r) + logger, t) (5.44) 

and 

loggKO,^) < logg+(0,r) +logg,^(r,t) (5.45) 

for all y & D, and all r,t > 0. For fixed k, let r = (1 — K,)t. Then plugging 
into ()5.44|) we have 

^loggf(0,t) > ^lz^logg,-(0, (1 - K)t) + ^loggK(l - ^)t,t). (5.46) 

We have already shown that 

lim ilogg7(0,t) = lim ylogg/(0,t) = -Ss{c). 

t— >00 J £-^00 t 

It follows from that 
- Ss{c) = limsupyloggf(0,t) 

- J™ n r:logg7(0, (l-K)t) + limsup — loggf((l-fi;)t,t) 

= —(1 — k)Ss{c) + Klimsup — logg|((l — K,)t,t). 

f— >oo K,t 

Thus 

limsup^loggf((l - K)t,t) < -Ssic). (5.47) 
An analogous argument with ()5.45p shows that 

liminf — logg^((l - K)t,t) > -Ss(c). (5.48) 
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Therefore, for each y G D, 

hm llogg|((l - /€)t,t) = -Ss{c) (5.49) 

almost surely with respect to Q. Now we need to show that for each uj, the 
limit is uniform in y G -D. By Lemma l5.2t there are constants Ki, K2, and 
K3 such that 



g+((l - ^)t,t) < - ^)t,t)+K,e 



It follows from ()5.49|1 that for any i/q fixed and any e > and 6' = 6/{l + e), 
we have 

liminf — logg7((l - K)t,t) > liminf — log (Kiq+{{1 - K)t,t) - K2e-^'''^''''' 

> liminf — log {Kiqf{{l - K)t, t) - 7^26"^^^''^'"'*' 

t^oo Hit \ 

> -Ssic) = -Ss/a+^c). (5.50) 
Since Ss{c) y S{c), inequality (j5.5(J|) implies that for any 6 > and h > 0, 



M^Piv'M)eUsic))>e 



-K.t(S{c) + h) 



for t sufficiently large. This proves part (iii) of Definition 14. II for k G (0, 1]. 
Using Lemma 1^721 and (|5.49|) we also have 

limsup^logg+((l - K)t,t) < hmsup ^log (K,q7^^.,{{l - K)t,t) + K^e-""''"'"^"'") 

< hmsup i-log (K,qy^ {{1 - K%t) + K,e-^^^'''^''') 

t— >oo ^ ' 

< -S{^Y+e)5{c). (5.5L 

Then letting 5^0, S(i+,)5{c) / S{c), so that ^fA^ and imply 

P{r]\{Kt) G Us{c)) < e-'^*(^(^)-'^) 

for t sufficiently large and 6 sufficiently small. Now using Lemma 13.21 as in 
the case of k = 1, part (ii) of Definit ion 14 . 1 1 follows . This completes the proof 
of Proposition 14.11 □ 



26 



6 Computing with Variational Formula 



In this section, we use formula ()1.7p to compute the propagation speed c*. 
We also derive some analytical bounds on c* and compare them with our 
numerical results. In our numerical simulations we consider a specific case 
of b{y,t) = bi{y)b2(t), y G R^, where bi{y) is a smooth periodic function of y 
and b2{t) is stationary Gaussian field. For 62 we use the Ornstein-Uhlenbeck 
process which solves the Ito equation 

db{t) = -ab{t)dt+ rdW{t), t>0, (6.1) 

where W{y,u!) is the standard Wiener process, X{0,u) = Xq^u) is a Gaus- 
sian random variable with mean zero, and variance p = r^/(2a). This is a 
mean zero stationary Gaussian process with covariance function 

EQMtMs)]=T{\t-s\) = ^e~'^\'-^\. 

In most computations, we will choose r = ^j2o?l^ so that the covariance is 

i^Q[fe2W&2(s)] = v^e-"l*-^l = r(|t-s|) . (6.2) 

By this choice of r, the l? norm of Vi^z) remains constant as a changes, 
so that the total energy in the power spectrum of the 62 remains constant. 
Because b\{%j) is periodic and because of the rapid decay of the covariance 
function V ^ it is easy to verify that b{y,t) = bi{y)b2(t) satisfies the assump- 
tions A1-A5. Assumption A6 follows from the Lipschitz continuity of bi{y). 

6.1 Numerical Computation of /i(A) 

To compute /i(A) = A^/2 + /'(O) + p(A) fsee 12. 2|) . we discretize the auxiliary 
initial value problem 

0i = ^A</.-A6i(y)62(t)0, 

0(2/, 0)^1 (6.3) 
using the Crank-Nicholson scheme: 



At 2' ^"^ -rm. . 2' 
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where D'^cf)^ denotes the standard second order discretization of the Lapla- 
cian (or \^y) centered at the discrete point {ym,tn)- The term F,J^ corre- 
sponds to the reaction term \hi{y)h2{t) evaluated at discrete points {ym,tn)- 
This scheme is imphcit, second-order in both time and space. In all simula- 
tions we use bi{y) = 6sm{67iy), where 5 > is a scaling parameter, and we 
compute on the domain y G [0, 1] with discrete grid spacing Ay = 0.01. We 
use an adaptive time step, since the implicit treatment of the reaction term 
requires that At < 2/{F^~^^). To generate realizations of &2(^) we integrate 
the Ito equation (j6.1|) using an implicit order 2.0 strong Taylor scheme (see 
fn\ ) with a discrete spacing Atb^ < 0.1 (At), where At is the adapting time 
step for the PDE evolution. 

To approximate /i(A), we iterate (j6.4|) for a very long time t = Tj and 
approximate 

/i(A) ^ /x,(A) = AV2 + /(O) + i log(||0||i). (6.5) 

By the results of the previous sections, fit converges to fi almost surely, so 
we need only generate one realization of the process. Alternatively, we could 
generate realizations of 62(^)5 evolve ()6.4|1 . and approximate 

MA)^i?Q[/.t(A,cI;)] = AV2 + /'(O) + EQ[^log(||0||O] 

1 ^ 1 

= AV2 + /'(0) + -^-log(||</.,||i). (6.6) 

In practice, we observe that Ty can be chosen much smaller when using ()6.6p 
instead of (j6.5p . since the mean converges must faster than an individual 
sample. In Figure^ we show one realization of the approximation /ij(A) com- 
pared with the ensemble mean ()6.6p . After a very short time, the mean shows 
relatively little fluctuation compared to the individual realization /ij(A,(2'). 
The variance of /Xt(A), shown in Figure El decays like 0(l/t). However, the 
need to compute a large number of realizations in ()6.6j) makes this approach 
no less computationally expensive than evolving only one sample for a very 
large time. So, we generally use ()6.5|1 to approximate /i(A). Typically we use 
a final time of Tf = 30,000. Figure El shows the convergence of /x computed 
by ()6.5p . In Figure El we show the distribution of /ij at different points in 
time. For this simulation we used A^ = 40, 000 realizations. 
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Figure 1: One realization of pt = ^ log(||0(?/, t)||i) = fit - jl - /'(O). The 
nearly flat curve shows the sample mean, iV = 40, 000 realizations. 

6.2 Computation of c* 

Computing the speed c* is quite simple: we use the above method to evaluate 
/i(A) at different points on the curve A t— > and then minimize in A. Direct 
simulation of c* would require evolution of (jl.lj) . a nonlinear, time- dependent 
PDE in two space dimensions. Finding c* by computing the curve /i(A)/A, 
however, reduces the PDE computation to one dimension. The trade-off is 
that we now have a minimization problem in A, but this is easily accomplished 
with a standard algorithm [7j. From the representation ()10.2|) . one can see 
that the curve yu(A)/A has the following properties: 

Lemma 6.1 The infimum of the curve over (0, oo) is achieved at a 
unique point A* G (0, Aq] where Aq = \/ /'(0)/2. Moreover, there are no other 
local minima. 

Proof: This follows from the fact that fi{X) = + /'(O) + p(A) with p 
being convex in A and p(0) = (see discussion leading to ()10.15|) ). The point 
Ao is the value of A where the infimum of the curve A/2 + /'(O) /A is attained. 
□ 

Next we consider the scaling b{y, t) ^-^ Sb{y, t) and the resulting enhance- 
ment of the corresponding speed c* = c*{6). It is known j2H] that if b{y,t) 
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Figure 2: Variance of AV2+/'(0) + i log(||0(?/, t)||i), = 40, 000 realizations. 

is periodic in both space and time that c*{6) = c*(0) + 0{6^) for 6 small 
and c*{6) = c*(0) + 0{6) for 5 as 5 — > oo. The following proposition gives 
analytical upper bounds consistent with this asymptotic behavior. 

Proposition 6.1 (Bounds on c*) For all S > 0, c*{S) satisfies the bounds 

(i) c*{5) > c*(0). 

(ii) c*{6)<c*{0) + 6\\h\\MM- 
(in) c*{6) < c*{0)^l + 6^pi. 

(iv) c*(5)=c*(0) if h{y,t) = h{t). 
From (Hi), we also have 

c*(5)<c*(0)(l + ^)+O(5^) 

when 6 is small. 

Proof: The first bound follows from ()10.15|) and the formula 

^ ^ A>o A - A>o 2 A ^ ^ 
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2.5 3 
xio" 



The function ip = log(0) satisfies 

i^t = ^Aij + ^\Vy^\'-\biy,t) 

Integrating ()(j.7j) over D x [0,t], we have 

^ / ^{y,t)dy = ^ [ [ \Vy^\Uydt-j 
Now let if: oo: 



JD 



p(A) = hm i / ip{y,t)dy > -A hm ^ / / b{y,t)dydt 



-XEr 



b{y,t) dy 



D 



0, 



(6.7) 



b{y,t)dydt. (6.8) 



almost surely with respect to Q. If b{y,t) = b(t), then the first integral on 
the right hand side of ()6.8p vanishes since iVj;'?/'^ = 0. Then taking the limit 
as t — s> oo we have equality: 



P(A) = Eq 
Hence c*{6) = c*(0). This proves (iv 



Sb{t) dy 



D 



0. 
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Figure 4: Distribution of j log(||0(y, t) ||i) at different times. 



For the linear upper bound (ii), note that 

llog^[gA<5/j6i(W(s))62(t-^)<is| ^ llog^[gA5||bi||^/o'|62(t-s)|rf^] 



\X\6\\h\ 



\b2{s)\ ds. 



As t — > cxo, this last term converges almost surely to |A|5||6i||£'q[|62|]- There- 
fore, c*{5) always satisfies the linear upper bound 



c*(.)^mf/f) < mf^.m,„„,„^,[„,, 
= c*{^) + 6\M\EQ[\h,\]. 



(6.9) 



Finally, for the quadratic upper bound, observe that under the scaling h i— *• 
the constant pi defined in A5 can be replaced by pi ^ A^5^pi. Then by 

(inrfii and (unm, 

p(A) < V2X^5^pi 
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and 



c*{5) = inf 



/i(A) 



A>0 A 



< 



. , A /'(O) X^6^pi 

mf - + ^^-^ H i-^ 

A>o 2 A 2 

2v/(l + ^2pi)f(0 )72 
c*(0)v/(l + 5V) 



c*(0)(l + 



c*(0) 



(6.10) 



□ 



If we make some additional restriction on the form of b{y, t) we also have 
a linear lower bound on the growth of c*{5) as 5 oo. 

Proposition 6.2 (Linear growth of c*) Let b{y,t) have the form 

N 

b{y,t) = J2b{iy)m) 
i=i 

where h\{y) are Lipschitz continuous and periodic in y, and bKt) are station- 
ary centered Gaussian fields such that the Assumptions A1-A6 are satisfied 
for b{y,t). Then the constant C G [0, +oo) defined by 

liminf^-^ = C (6.11) 

S—fOc 

is equal to zero if and only ifb{y,t) = b{t). 

Proof: The fact that C G [0, +oo) follows from Proposition 16.11 Also, if 
b{y,t) = b(t) then (5 = since c*{5) = c*(0) for all 5 > 0. By Lemma O 
there is a unique A = A^ G (0, Aq] such that 

c*(5) = inf^ = 4^. 
A>o A A^ 

Let Sj — >• oo as J — > oo and suppose that \imsup j^^{Xs^6j) < M. This 
implies that 
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So, in this case the result holds with C = /'(0)/M. 

Now suppose Xs^Sj is unbounded as j — oo. By Proposition llU.ll there 
is a positive constant K such that 

piXs^5j) > KXs^5j > 

for j sufficiently large. Note that Proposition 110. ll treats the case of 5 = 1; 
this is why we use p{XsjSj) instead of p{Xsj). Therefore, 

hmmf — 7^ = limmf — ^ + ^^-^ + > K > 

j-rco dj j^oo 2oj XsjOj Xs.Oj 

since A^^^ G (0, A*] and Xs.6j oo. Hence C > K > 0. □ 

The results of our numerical computations suggest that c*{5) is a mono- 
tone increasing function of 6, and they conffim both the quadratic and linear 
growth of c*{6) for small and large 6, respectively. In Figure El we plot 
c*{S) — c*(0) on a log- log scale for a few values of the covariance parameter 
a. We observe a transition form quadratic scaling in the small S regime to 
linear scaling in the large S regime. We also plot the upper bounds gi and 
(?2 given by (ii) and (iii) of Proposition 16 . II using a = 16.0. We observe that 
these bounds are at least an order of magnitude greater than the numerically 
computed enhancement. 

Next, by varying the parameter a in ()6.2j) . we consider the effect of tem- 
poral correlations on the enhancement of the speed c*. Because of the simple 
structure of the field bi{y)b2(t) we can compute the constant pi appearing in 
Assumption A5 and Proposition 16.11 Using (j6.2p . 

supr(i/i,i/2,0,r) = IMlEgMOMr)] = v/^||6i||^e-"^ 
yi,y2 

Hence pi = Therefore, from (iii) of Proposition 16. H 

c* < c*{0)^{l + 6^a-y^h\\l) (6.12) 

So, as a ^ oo, c* c*(0). This limit corresponds to the correlation length 
1/a becoming very small and is consistent with the case of periodic time 
dependence: faster temporal oscillation of the shear tends to decrease the 



34 



10' 



10' 



.10 



10" 



10 



■ o 


a=1.0 


■ ■□ 


a=4.0 


■ + 


a=8.0 


A 


a=16.0 


V 


g^, a=16.0 




Qj, ce= 16.0 



10" 





10 



10 



10 



10 



10 



Figure 5: Log-log plot of front speed c*{5) versus S. For comparison, the 
solid lines on the bottom have slope p = 2.0 and p = 1.0. The data sets gi 
and g2 represent the upper bounds ()6.9|) and ()6.10|) . respectively, contained 
in Proposition 16.11 Both bounds were computed for a = 16.0. 



enhancement of the front speed [201 12Z1 12H|- As a ^ 0, the bound ()6.12p 
blows up. From (ii) of Proposition 16.11 however, we also have 



c* < c*{0) + 6\\h\\EQ[\b2\] 
= c*{0) + 6a'/'\\h\\EQ[\Z\] 



(6.13) 



where Z is a normally distributed random variable. This implies that c* 
c*(0) as a — > 0, as well. So for a G (0, cxd), there must be some optimal 
correlation length 1/a so that enhancement is maximal. 

For a G [1/10,2^^/10], we computed the expected speed c*{6) for fixed 
amplitudes 6 = 2.0 and 5 = 40. Note that for each a, we must choose the ini- 
tial points 62(0) to have variance E[b2{0)'^] = ^Ja so that the process remains 
stationary for each a. Also, as a becomes large, we adjust the PDE time 
step so that At < 0.5/a, in addition to other restrictions already mentioned. 
Otherwise, the numerical method cannot resolve the fast oscillations of the 
shear process, and the speeds diverge as a grows. 

In Figure El we plot c* versus a on a log-log scale. We also plot the upper 
bounds g\ and g^ given by ()6.13|) and ()6.12|) for 5 = 40. In general, these 
bounds are rather coarse, an order of magnitude larger than the numerically 
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Figure 6: Log-log plot of c* vs. covariance parameter a {5 = 2.0,40.0). 
For comparison, the solid lines have slope p = 1/2. The data sets gi and 
g2 represent the upper bounds ()6.9|) and ()6.10|) . respectively, contained in 
Proposition 16. II Both bounds were computed for 5 = 40.0. 



computed enhancement, but as a — oo, the bound g2 lies relatively close 
to the data. As a — >• 0, gi is the better bound, predicting that c* c*(0), 
although the scaling of the bound is different from the scaling observed in the 
data. By computing the slope of a best-fit line through the points on the log- 
log plot, we find that for 6 = 2.0, c* scales according to c* = c*(0) + 0(a°-^°) 
and for S = 40, c* = c*(0) + 0(a°-^^). So in both cases, \c* - c*(0)| 
much faster than the 0(a°'^^) convergence predicted by the bound gi in 
fl6.13|) . This scaling behavior can be understood by analogy with the periodic 
case. Note that 5bi{y)b2{t) = a^^^5bi{y)b2{t) where 62 has unit variance 
and correlation length 1/a. If b{y,t) = bi{y)b2it) and 62 is periodic with 
very long wavelength, then the enhancement is approximately equal to the 
enhancement caused by bi only (a steady shear). The very slow oscillations in 
the shear field do not significantly slow the front. In the random case, when 
a is small the correlation length is very large. So by analogy, we expect that 
in the random case, when a is very small, c* will behave as if the shear were 
just a^^'^Sbi- In the small amplitude regime, c* scales quadratically with 
amplitude. Hence c* ^ c*(0) + 0{{a^/^6f) = c*(0) + 0{a^/^^), which is 
consistent with our numerical computations for a small. Figure [7| shows the 
data for 6 = 40.0 in terms of correlation length 1/a. 
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Figure 7: Dependence of c* on correlation length 1/a {5 — 40.0). 

7 Conclusions 

We have considered the propagation of KPP reaction fronts in temporally 
random shear flows inside an infinite cylinder. We showed that, under as- 
sumptions A1-A6 on a Gaussian shear field, the fronts speeds obey a varia- 
tional formula that extends the known variational formula in the case of pe- 
riodic media. We performed analysis and computation of front speeds based 
on the variational principle. When the shear field is Ornstcin-Uhlenbeck in 
time, we numerically demonstrated the quadratic and linear speed growth 
laws in the shear root mean square amplitude, and speed dependence on the 
shear temporal correlation length. We also derived basic bounds on the front 
speeds and compared them with the computed speeds. Developing methods 
to improve these bounds will be left as a future work. 
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9 Appendix A: Estimates on X^^(s) 

In this section we derive some technical estimates on the family of processes 
that follow from our assumptions on the field b and the Borell in- 
equality for Gaussian fields. In particular, we prove Lemma 13.21 which is 
needed for the bounds ()3.9j) and ()4.15|1 . Let use first note that by changing 
variables r = s — t, v = s + 1, it is easy to see from Assumption A5 that 

/ snp T{yi, 1/2, s,t) ds dt < 2 / T{r)drdv 

Jo yi,y2 Jo Jo 

< 2V2piT (9.1) 

and for H e [0,T], 

j j supT{yi,y2,s,t)dsdt < V2\T - H\ / t{r)dr 
Jo Jh yi,y2 Jo 

< V2\T-H\pi. (9.2) 

Let p{s) e C([0,+oo),R"-i) with p(0) = be fixed. For y e D, define 
Pyi^) = 2/ + P{^)- For fixed t > 0, the integral 



f{y,s)= / b{py{T),t-T)dT (9.3) 
Jo 

is a Gaussian random field over M = D x [0,t\, with respect to the mea- 
sure Q. The Borell inequality for Gaussian fields states that if ||/|| = 
sup(^ g-jgj^^ f{x, s) is almost surely finite, then < oo and for any u > 0, 

Qm\-mf\\]>u)<e'i (9.4) 

where = sup(^, (see Q). By dHIH), (x^ < 2V2pit. Using in- 

equality ()9.4j) . we can control deviations of ||/||, if we bound the growth of 
E\ 



Lemma 9.1 There is a finite constant C > such that 

E[\\f\\]<Ct'/'. (9.5) 
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Proof: The expectation can be bounded by the metric entropy rela- 

tion PP 

^[ll/ll]<C [\\ogN{e)f' de 
Jo 

where 6 = diam{M) /2 in the metric 

di{x,s),{y,z)) = E[if{x,s)-f{y,z)YY^' 

and A^(e) is the minimum number of e balls required to cover M. Using ()9.1|) . 
f|9.2p . and Assumption A6, a straightforward computation shows that 

E[{f{x,s)-f{y,s)f] <C\x-y\t 

and 

E[ifiy,s)-fiy,z)f]<C\s-z\ 
for some finite constant C, independent of p. Therefore, 

d{{x, s), {y, z)) <C,{\s- z\f'^ + C2{\x - y\t)"' , 

and there is a constant C3 independent of t and e such that d{{x, s), {y, z)) < e 
whenever |s — 2;| < C^e'^ and — ?/| < For e G {0, diam{M) /2], we have 
the bound ^ 

iV(e) <max(C5^,l) 

and 

^[||/ll] < C I {\og{C,-^)y/' de 

Jo ^ 

= C,f/' f\\og[\)f/\de<C,e'\ (9.6) 
Jo £ 

□ 

Note that the constants depend on the assumed properties of the process 
h and the size of the domain D, but not on the particular function p{s). If 
u > 2C7t^/^, then by ^U^ . 

Q (11/11 >U)< e-(-^ll/ll)V2<x? < g-«V8<x? < g-«V8pit_ 

It now follows that 
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Lemma 9.2 For any rj > and for t > to — to{r]) — {2C-j/rff', 

Q I sup f b(y + p(r) ,t-r)dr>r]t] < e'^^^/^P' . 
\yeD, se[o,t] Jo J 

for any p e C([0, oo), i?"), p(0) = 0. 

In applying this lemma, the continuous function p will be a realization of the 
Wiener process W2^(s). 

Lemma 9.3 For r] > 0, z e R^, define the Markov time 

r^^,{t) = mm{s > 0| \X^'\s) - x\ > r]t}. 

with T^^zif) = +00 if the set on the right is empty. Then there are constants 
Ki, K2 such that 

Q [p ( mj^r,,,(t) <t^> e-^^'''*/^^ < K^e'''^^""'' 

for all t > 0. 

Proof: Note that for the Wiener process Wi{s) with IVi(O) = 

P I sup >r]t\ < 2J- f e-^'/2 

< K^e-""'/\ (9.7) 

The point of the lemma is that at large times, and almost surely with respect 
to Q, the process X^'*(s) still behaves like a Wiener process, even though the 
drift term b{y,t) is not uniformly bounded in t. By definition of T^,z(t), 



P ( inf T,,,{t) <t] = pI sup Iftiv, s) + W: 



i{s)\>Vt\. 
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Using Tchebyshev's inequality, (j9.7p . and Lemma we see that for any 
T] > 0, a > 0: 

q(p( M Tr,,,{t)<t]>a] < a-'EQpl sup \ftiy, s) + W^is)\ > rit 
\ V^e^" / / \se[o,t],2ei?" 

a-'EpQ ( sup \ft{y, s) + Wi{s)\ > r]t 

< a-'EpQ{ sup \ftiy,s)\>r]t/2 

se[o,t],zei?" 

+a-^P[ sup \Wi{s)\>7]t/2 

ys(^[0,t],zeR" 

for t sufficiently large, for some constants Ki, K2 > 0. The result now follows 
from a choice of a = e'^^"^ ^1'^ . □ 



Corollary 9.1 There are constants Ki,K2 > such that, except on a set of 
Q -measure zero, 

sup P (r^,,(t) <t)< Kie-^''''' (9.8) 
for t sufficiently large depending on Co and rj. 

Proof: Note that this is Lemma f3. 21 in the case that k = 1. Lemma f9. 31 and 
the Borel-Cantelli lemma imply that outside a set of Q-measure zero 

P inf^ r^,,(A;) < A;^ < g-^^''''^/^ ^g g) 

if G Z is sufficiently large. Now we want to extend this to all real t 
sufficiently large. Let t E [k, k + 1], t = k + r, r G [0, 1]. 

sup P sup |X^'*(s) - xo| > tr] 
zeR^ yse[o,t] 

< sup P f sup |X^'*(s) - xol > t7]/2 I + 

+ sup P I sup |X"'*(s) - X''*(r)| > tri/2 

zeR" \se[T,t] 
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By the Markov property, this is bounded by 

< sup P I sup \X''\s) -Xo\> t7]/2 ] + 
zeR" Yse[o,T] 

+ sup P I sup |X^'^(s) - xo\ > tr]/2 | < 

< P f sup \X''\s) - xo\ > kri/2 ] + 

yzeR" ,te[k,k+l],selO,l] 

+P I sup \X''\s) - xo\ > kri/2 1 (9.10) 
Yfei?",se[o,fc] J 

By ()9.9|) . the second term on the right side of ()9.10p is bounded (Q-a.s.) by 

P I sup \X''\s) - a;o| > A;r//2 | < 6"^^''''= (9.11) 

for k E'L sufficiently large. To bound the other term in ()9.10p . it suffices to 
show that 

P I sup \fk{y.s,r)\ > kr]/2 ] < e'^^^'^^ (9.12) 

Yy6D,re[0,l],se[0,l] 

for k e'L sufficiently large, where 



fk{y,s,r) 



I h{W^{T),r + k-T)dT. 
Jo 



Note that fkiy, s,r) is a centered Gaussian field over D x [0, 1] x [0, 1], and 
its distribution is invariant with respect to k > 0, due to the stationarity of 
b{y,t). For any fixed path iyj^(a;), the Borell inequality implies that for k 
sufficiently large 

Q ( sup \fk{y,s,r)\ > kr^/2] < K,e-^^^''' 
ys/GD,re[o,i],se[o,i] J 

for some constants K^, Kq > 0, independent of k and the realization (u;). 
Therefore, proceeding as in the proof of Lemma 19.21 we see that 

Q 1 P( sup \fk{y,s,r)\ > kv/2) > e-""'^"'"/^] < K^e-'''^^"^''\ 

J/G-D,re[0,l],se[0,l] 
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Now (|9.12|) follows from the Borel-Cantelli lemma. We complete the proof 
by combining ^TH^ and (jHUS)- □ 

Proof of Lemma 13. 2t We have just proved Lemma 13.21 in the special 
case that k, = 1. For k < 1, modify the preceding bounds for the field 



f{y,s) = / b{py{T),t-T)dT 

Jo 

considered over = D x [0,/tt]. Now we have erf = sup(^^^g-^^j^,j Eq[P] < 
Pint, so we find that < C^/kI for some constant C > 0. Then, just 

as in Lemma 1^731 we have 

Q\P\ sup \X'\s) -x\>^t\> e-^2ry2*/2- ) < Kie-^^'''*/^^. 

The rest follows as in the preceding corollary. □ 



10 Appendix B: Lyapunov Exponent /i(A) 

In this section we prove Lemma im assuming that A1-A6 hold for the random 
field b{y,t). We study the limit 



fi{X,z) = /'(0)+limllogE 

= /'(0) + ^ + limilogE, 

= /(0) + y+p(A,l/). 
By the Feynman-Kac formula 



Xj;b{W^is),t-s)ds 



(10.1) 



p{\y) = lim -log0(y,t) 

t— >oo t 

where > solves that auxiliary initial value problem 



(10.2) 



<f)t = -A(f)-Xb{y,t)(j) 



(10.3) 
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The equation (|l(J.3p is called the parabolic Anderson problem (see (TU] and 
[12]). The proof that p(A) exists (and /i(A)) almost surely with respect to Q, 
independent of relies on the subadditive ergodic theorem and a recent result 
in il2|, provided we assume the necessary decay of time correlation of the 
process b{y,t). For simplicity of notation we ignore the dependence on —A, 
—Xb (-^ b. Following [T2], we define for any continuous path X G C{[0,t],R) 
the exponential 

For any fixed path X, C,{t,X) is lognormal with mean and variance 



EQm,X)] = e-, Var[^{t,X)] = e'^ - ^ 



where 



I [ T{Xs,Xr,s,r)dsdr <2V2pit, (10.4) 
Jo Jo 



by (jy.ip . Note that is bounded independently of the particular path X. 
For < s < t, define the random variables 



qi{s,t) 
(ls{s,t) 



sup 



Ey[e^O-'biW(r),t-r)dr^_ 



Using the subadditive ergodic theorem, we will show that the limits 



and 



lim -logg/(0,t) = pi 



lim \ \ogqs{^,t) = Ps 



;io.5) 



:io.6) 



exists and are finite, almost surely with respect to Q. Then we will show 
Pj = Ps- By the Markov property of the Wiener process we have for and 
s < z <t: 



qi{s,t) 



miE. 

y 

iniE. 

y 



> infK, 



g/o'-^ fe(M/(r),t-r) dTgj/r; b{W{T),t-T) dr 
elt' b{W{r),t-r) dr-^^^J^Z; biWir),t~r) d^|yf/^_J 
e/o b{W{r),t-r) dr ^^^^-^ (r'),.-r') dr' 



qi{s,z)qi{z,t). 
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Therefore, log(g/(s,t)) is super-additive: 

log(g/(s, t)) > log(g/(s, z)) + log{qi{z, t)) 

for any < s < z < t. Similarly, the function \og{qs{s,t)) is subadditive. By 
the stationarity of b, 

Tz log(g/(s, t)) = log(g/(s + z,t + z)) 

for any z > 0. Moveover, log(g/(s,t)) is integrable: 

EQ[log{qi{s,t))]<^ogEQ[qj{s,t)] < loginf EpEQ[e^r '>mrU-r)dr^ 



and 



< loge'2 < V2pi\t — s\ 



EQ[logiqr{s,t))] > Eq llogiEpe'^'^ ^(^Hr),t-r),r^ 



(10.7) 



> EpEq 



inf / b{Wy{T),t-T)dT) 

.y Jo 



(10.8) 



This last term is finite, by the Borell inequality. Note also that 

sup |logg/(s,t)| < sup \b{y,t)\, 
s,te[o,i] te[o,i],yeZ) 

and the latter is integrable with respect to Q. It now follows from the sub- 
additive ergodic theorem (Theorem 2.5 of j2I) and the continuity of q{0,t) 
with respect to t that the limit ()10.5|1 exists almost surely and is finite: 

lim 1 loginf t) = lim l^SMM) = ,„p E,V°^MOM ^ 

i^oo t y t^oo t t t 

Also, by (fTITfl) . Pi < v^pi. 

In order to apply the ergodic theorem to (1/t) log(g5(0, t), we need to 
show that (1/t) log(gs'(0, t) is integrable. By Jensen's inequality. 



EQ[\ogqs{s,t)]>snpEQEy[l b{W{r),t - r) dr] = 0. (10.10) 

y 



45 



Also, it follows from the Borell inequality and Theorem 3.2 of (p. 63, let 
a = 1) that there is a finite constant Kq > such that 

EQe''''PyIo~^biWy{r),t~r)dr ^ < OO (10.11) 

if (T^ < i . Thus by (llU.4p . there is a constant Ki > such that (|1U.11|) 
holds when |t — s| < Ki. Now for any s < t, let be the smallest integer 
greater than \t — s\/Ki and s = to < < ^2 < ■ ■ ■ < ^Af = ^ with — tj| = 
At = |t - s|/A^ < Ki for alH = 0, . . . , A^ - 1. Using the subadditivity of 
\og{qs{s,t)) and Jensen's inequality, we see that 

EQ\og{qsis,t)) < ^EQlog(g5(ti,ti+i)) 

i=0 
N-l 



< ^log(EpEQe^ 



i=0 

By ()10.11|) . this right side is bounded by A^log(A'o) < oo. Moreover, 

sup \\og qs{s, t)\ < K, sup \b{y,t)\ (10.12) 

s,t<^[0,Ki] t<^[0,Ki],y<^D 

which is integrable. Hence, we can apply the subadditive ergodic theorem to 
conclude that the limit 

t^oo t t t 

holds almost surely with ps a constant, ps G [0, oo). Note that convergence 
along continuous time follows from ()10.12|) . the continuity of ^^(0,^:) and 
Theorem 2.5 of 

Clearly p/ < ps- To show that pi = ps, we will need a kind of Harnack 
inequality to compare the quantities qi{0,t) and qs{0,t). Such a result has 
been obtained in in the case that b{y, t) is Gaussian in both space and 
time, with a white-noise temporal dependence. Under the assumptions A5- 
A6, however, the arguments of |12| imply that the following estimate also 
holds in the present case. 



Theorem 10.1 il2^ For any fixed M > 0, there are positive constants ci, C2 
such that outside an event of Q— probability e"*"" , one has 

inf Ey [e(n,X)] > cie-^^""''" fsupE, [e(n,X)] - 6"^"''') . 

\y&D J 
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From this result it follows immediately that 

lim - log inf (f){y, n) = pj = ps, 

n— >oo n y 

and by ()10.9p and ()10.13|) . we see that this extends to continuous time 

hm iloginf = pj = ps = piX) - /'(O) - ^. (10.14) 

i^oo t y Z 

We have now shown that p(A) (and thus /i(A)) in ()10.1|) is well-defined, 
independent of z E -R", almost surely with respect to Q. 

The convexity of p{X) follows from the convexity of p(A). Let r G (0, 1), 
Ai, A2 G R. By Holder's inequality, 

Applying this to ()1().1|1 . we conclude that 

p(rAi + (1 - r)A2) < rp(Ai) + (1 - r)p(A2). 

Clearly p(0) = 0. Since b{y,t) has the same distribution as —b{y,t), we see 
that p(— A) = p(A) for all X E R. Hence p and p are even functions of A. 
Since p(0) = and p is convex and even, we conclude that 

p(A) > and p(A) > AV2 + /'(O), VA G R. (10.15) 

This completes the proof of Lemma 11.11 

We conclude this section by demonstrating that p(A) grows linearly with 
A as A 00. 

Proposition 10.1 Let b{y,t) have the form 

N 

b{y,t) = J2b{{y)m (10.16) 

i=i 

where b\{y) are Lipschitz continuous and periodic in y, and b^t) are station- 
ary centered Gaussian fields such that the Assumptions A1-A6 are satisfied 
for b{y,t). Then there is a constant K > such that for X sufficiently large, 
p(A) > KX. 
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Proof: In the case that b{y, t) is a Gaussian field with white-noise time 
dependence, the authors of [T2| studied the behavior of p(k) as k — > 0, where 
K > is a diffusion constant (replace A/2 with kA/2 in p0.3|) ). Here we 
modify some of their ideas to treat the large advection limit when b has the 
form mj^ . 

Let At be the set of continuous, piecewise linear functions g such that 
g{0) = 0, g is linear on the intervals [it/k, {i + and 

g{{i + l)t/k) - g{t/k) = ±t/k, i = 0,1,2, k - 1 

where > is a large integer that will depend on t. Using the subadditivity 
arguments of |T2j one can show that there is a constant Ci > such that 



Q f sup [ b{f{s),t-s)ds>{ci~e)t)>l-e (10.17) 
\feAt Jo J 

if k is chosen to be sufficiently large, depending on t. This implies that by 
choosing e small there is a set of probability at least 1/2 such that we can 
find /(s) = /(s,a>) G At satisfying 

/ b{f{s,u;),t-s,u)ds>^, (10.18) 
Jo ^ 

and we expect that Brownian paths staying close to / will make a significant 
contribution to the exponential in the definition of p(A). For a constant 7 > 
to be determined and / G At, we let Bt{f,'j) be the 7-neighborhood of / in 
C([0,t],i?"-i): 

5,(/,7) = {X G C([0,t],i?"-i) I ||X - /lico < 7} . 

Then there are constants Ki, K2 independent of t and of / G such that 

P{Bt{f,j))>K,e-^^('-^'/^'^' 

for t > 1 (see J2])- Now using assumption A6, we see that for any path 

XeBtif,^), 

ft ft ^ ft 

1/ b{X{s),t- s)ds- b{f{s),t-s)ds\ <-fMj2 W2is)\ds (10.19) 
Jo Jo Jo 
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where M is the maximum of the Lipschitz constants for the functions {b{{y)}j^i 
By (I1U.18|) and (jl(J.19p with e > sufficiently small, there is a set of Q- 
probability at least 1/2 such that 



,A f*b{W^(s),t-s)ds 



> Er 



,A j'^b(WUs),t-s)ds 



where V = Io\^2i^) \ / G A is chosen to satisfy (|1(J.18|) . For t 

large, independently of A and 7, V can be bounded by 



N 



v<t(^E[\hm] + i) 



except on a set of probability less than 1/4. Therefore, we can choose 7 small 
so that 

7< 



Cl 



4M(ELi?[|&K0|] + l) 



Then by choosing A large we obtain from ()10.20j) 



E, 



,XSl^h{WUs),t-s)ds 



> e 



Xcit/8 



with Q-probability at least 1/4, for t sufficiently large, independently of 
A. Since the limit defining p(A) exists Q-almost surely, this establishes the 
lemma with K = ci/8. □ 
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